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Abstract 

We calculate the one-loop effective potential at finite temperature for a system of 
massless scalar fields with quartic interaction A0 4 in the framework of the boundary 
effective theory (BET) formalism. The calculation relies on the solution of the classical 
equation of motion for the field, and Gaussian fluctuations around it. Our result is non- 
perturbative and differs from the standard one-loop effective potential for field values 
larger than T / VX 



1 Introduction 

In thermal systems containing massless bosons, a direct implementation of perturbation 
theory for the calculation of thermodynamic quantities is problematic. Vanishing Matsubara 
modes bring, then, infrared divergences that render the naive perturbative series essentially 
non-convergent and meaningless [1]. Therefore, the only way to proceed with a sensible 
perturbative calculation in this realm is to reorganize the diagrammatic series by resumming 
certain classes of diagrams. In particular, in the early days of the computation of the one- 
loop effective potential, it was immediately realized that one should resum the so-called ring 
diagrams even in the limit of weak coupling, and that this recipe provides a thermal mass 
for the bosonic (thermal) propagator [2]. Resummations of this sort can be performed and 
improved in different ways. We refer the reader to the reviews of Refs. [3, 4, 5] and to Ref. 
[6] for a longer discussion and a list of more specific references. 

In a recent paper [7], we have proposed an alternative approach to thermal field 
theories, which we denoted boundary effective theory (BET). In the present paper, we 
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apply the BET formalism to compute the one-loop effective potential at finite temperature 
for a system of massless scalar fields with quartic interaction A</> 4 /4!. The calculation relies 
on the solution of the classical equation of motion for the held, and Gaussian fluctuations 
around it. Our result is non-perturbative and differs from the standard one-loop effective 
potential [2] for field values larger than T j\f\. 

The boundary effective theory (BET) is a natural way to organize the calculation 
of the partition function of a quantum system at finite temperature, where one slices the 
functional space of /3-periodic fields into sectors where the boundary value of the field is 
fixed [7] . More precisely, this amounts to writing the partition function Z as 

Z = J [D<j> {xj\ pp[Mx),M^)] > (!) 

where 4>o(x) is the field at the boundaries of the imaginary time interval, and pp[<fio(x), (f>o(x)] 
is the functional density matrix diagonal element, given by the functional integration over 
all the fields <f>(r, x) that have this particular value at the time boundary, 

P/3 [<M*),<M*)] = J [D4>{t,x)\ e- s *M, (2) 

<t>(0,x)=4>{P,x)=4> (x) 

where S E [(f)] is the Euclidean classical action. 

The functional density matrix diagonal element pp[<fro(x), (f> (x)} is an expression in- 
volving only the static boundary field 4>o( x ), which already contains all the temperature 
dependence. Thus, correlations calculated within this reduced theory encode all information 
about the thermal distribution of the fields <fio( x )- One can define a dimensionally- reduced 
action Sd[4>o(x)] through 

Pfi [Mx),M*)]=e- s < [M " )] ■ (3) 



The reduced theory also contains the relevant infrared physics. Indeed, the double 
integral structure of the partition function naturally separates the static modes <fio(x). As a 
consequence, the reduced theory amounts to a resummation of an infinite class of diagrams 
of naive perturbation theory. This was verified in Rcf. [8], where the effective action 
obtained in [7] was used to compute the pressure to lowest order in the BET formalism. The 
resummation of ring-like diagrams emerges directly from the theory, and the corresponding 
pressure is in good agreement with recent calculations using weak-coupling and screening 
perturbation theory [9, 10, 11]. 

One should notice that BET is not a particular case of a procedure known in the 
literature as Dimensional Reduction (DR) (see Refs. [12, 13, 14, 15, 16, 17]). Indeed, 
DR methods produce effective theories for the zero Matsubara mode and, as such, are 
high-temperature approximations in character. The procedure that we utilize (BET) yields 
an alternative dimcnsionally-reduced effective theory for the physical field (j>o(x), and it is 
essentially different from DR. 

The aforementioned calculations using the BET formalism are non-perturbative, the 
results being expressed in terms of the classical solutions 4> c of the Euler-Lagrange equation. 
When these classical solutions are computed exactly, instead of being expanded perturba- 
tively, they automatically resum the infinite series of tree diagrams in the strong field regime. 
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Note that the Euler-Lagrange equations must be solved with non-trivial boundary conditions 
at the endpoints of the imaginary time interval: </> c (t, x) must be equal to <j>o(x) at the time 
boundary. To emphasize this functional dependence, we will denote the classical solution 
by <j) c [(f>o(x)}. In Ref. [7], it is shown that the effective action for the fields 4>o(x) admits a 
simple expression in terms of the classical solution (j) c [(j) (x)]. In order to extract the explicit 
dependence on (f>o(x), one must solve the Euler-Lagrange equation for arbitrary boundary 
conditions, which is not a feasible task in an interacting theory (but can, in principle, be 
done numerically). 

In this paper, we solve the full classical equation with constant boundary conditions 
0o (x) = 4>o, allowing for the calculation of the one-loop effective potential of the theory. The 
complicated non-linear dependence of <fi c on the boundary field <fio leads to a non-perturbative 
result that should be a good approximation even at large fa. 

Finally, it should be clear that the use of the word "boundary" has no relation to 
other uses, often in a topological sense, such as in holographic gauge-gravity duality, etc. 

The structure of the paper is as follows: in Section 2, we revisit the usual prescription 
to obtain the effective action and the effective potential in the functional integral formalism; 
using the results of Ref. [6], we write the effective potential in terms of the solution of 
a certain differential equation; in Section 3, we apply the method to the massless A^ 4 /4! 
theory, and discuss the results for the effective potential; finally, in Section 4, we present 
our conclusions. 



2 The one-loop effective potential in BET 

Let us consider the following classical Euclidean action, 

1 



(4) 



where (d 4 x) E is a shorthand for drd 3 x. We assume that U(4>) is some single- well interaction 
potential. Following the standard procedure for obtaining the effective action, we couple the 
boundary field to an external current j(x), and define the generating functional for the 
reduced theory [18]: 

Z[j(x)} = J [DM*)} e- s 'l*°(«>] , (5) 

where 

Sj[fa(x)] = S d [fa(x)] ~ P J d 3 xj(x)fa(x) (6) 

is the action for boundary fields in the presence of j(x). The free energy functional can be 
obtained from the generating functional Z[j(x)] as 

F[j{x)} = -\ lim lnZ[j(x)}. (7) 

By performing a Legendre transform of F[j{x)] 1 one formally obtains the effective action 

r[(Mx))j] = F\j{x)] + J d? x j(x)(Mx)h ■ ( 8 ) 
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Its argument is the expectation value of the field in the presence of the external current. 
The index j in (4> (x))j is to stress the dependence of the expectation value on the external 
current. One can think of T as minus the pressure of the system in response to an external 
current j(x). It can be shown that, at one-loop order, we have F[j(x)} — T[((j) (x)}j], and 
that {(f>o(x))j is the saddle-point of Sj[<j>o(x)]. 

In Ref. [7], it was shown that the one- loop effective action for 4>o(x) is given by 



priM*)] = smm*)]] + ^Tr In (A^ 1 + U"(MMx)})) 



(9) 



where is the time-ordered thermal propagator and 4> c [4>o{x)](t, x) is the classical solution 
of the Euler-Lagrange equation for fixed time boundary value 4>o(x), i.e. 



n E ct> c (T,x) + U'(4> c (r,x)) = 0, 

<t>c(0,x) = <j> c (P,x) = (j) (x) , 



(10) 



with \3 E = — (<9 2 + V 2 ) the Euclidean D'Alembertian operator. We see that, when written 
in terms of <f> c [<f>o(x)], T[(f> a (x)] is the same functional as the one- loop effective action at zero- 
temperature. In terms of graphs, cq. (9) includes the diagrams represented in the figure 1. 
In this figure, the first two terms represent the classical action evaluated at the classical field 






Figure 1: Diagrammatic representation of the content of cq. (9). The black dots at the 
endpoints of the trees represent the boundary field <fro( x ) at which the effective action is 
evaluated. 



configuration <p c . The tree structure of <p c when expressed in terms of (f>o{x) is manifest in 
the Green's formula that relates (j) c to its boundary value <fro, 

Mt,x)+J dT'd 3 yG(T 7 x;T' 7 y)U'(Mr\y))^ J d 3 y Mv) [d T 'G(r, x; r', y)] ^ , (11) 

where G(r, x; t', y) is the Green's function of the Euclidean D'Alembertian defined by 

□»G(r,aj;r',y) = 8(t - r')5(x - y) , G(t, x; 0, y) = G(t, x; 0, y) = . (12) 

The last term in the figure 1 is the diagrammatic content of the term ^Tr ln(- • • ) in eq. (9). 
Obviously, the non-linearities in <fi are important only at large 0o - when fa is small 
(compared to the temperature), all the trees in the figure 1 simplify into their lowest order 



4 



term in (j> (which is linear in O and amounts to solving eqs. (10) by neglecting the non-linear 
potential U'{<j> c ))- 

In order to obtain the explicit dependence of r[0o(aj)] on 4>q{x), we must solve Eq. 
(10) for arbitrary boundary conditions, which is in general unfeasible for boundary fields 
with an arbitrary x dependence. A simpler, yet very useful quantity to compute, is the 
effective potential, which is essentially the effective action evaluated for uniform boundary 
field configurations, 



V eS (<f> ) = (T[<Po] - T[0])/V (for constant cj> {x) = <j> ), 



(13) 



where V is the volume. In this case, the classical solution (j) c [<f>o\ depends only on r. Defini- 
tion (13) is such that V e s(0) — 0. As shown in Ref. [7], the quantity r[0]/V is the negative 
of the free pressure, i.e. (— 7r 2 T 4 /90). 

The calculation of the term ^Tr ln(- • • ) in Eq. (9) was done in Ref. [6] as an inter- 
mediate step to obtain the pressure in the context of a semiclassical approximation. There, 
it was shown that this quantity can be expressed, for constant <fio, in terms of solutions of 
the equation for small field perturbations propagating on top of the classical solution (f> c 



Tr In (A; 



u"{<j> c \M)) 



V 



d 3 k 

(2tt) 3 



In [2 ( V ((3,k 2 )-l)] 



where r\ (r; k 2 ) is the solution of 



[d 2 -k 2 -U"(UM(r))]v(r 1 k 2 ) = Q 



r,(0;fc 2 ) = l 



dr] 



(0;fe 2 ) = . 



(14) 

(15a) 
(15b) 



Therefore, we obtain the following (non-renormalized) expression for the effective potential, 



f3VV eS (<Po) = S B [<p c [cj>o}} + ^ 



d 3 k 



ln[2(r7(/?,fc 2 )-l)] - 0T[O] 



(16) 



where the integral over 3-momenta is regularized by the introduction of a cut-off A. 

In order to rcnormalizc (16), we add the standard one-loop counterterms C.T. (ob- 
tained by using the same cutoff regularization as the one used in Eq. (16)), and subtract 
the zero-point energy term (3k/ 2 from the classical action, obtaining 



/3VV eSR {4> ) = S E [ct> c [M}+ lim ^ 

A— >oo z 



d 3 k 



(2tt) 3 



In [2( V ((3,k 2 )-l)]-f3k 



C.T. 



where the counterterms read 



C.T. = V 



Ci 



dT<t> 2 c {r) + V 



Co 



with 



a = * [ A (k) 

Cl -2j {2^ Ak) 



16tt 2 



dT<£(T) , 



A 



(17) 
(18) 



(19) 
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and 

a^u \2 / ^ 1 

32tt 2 V"*M 2 

where /U is the renormalization scale. 



C ^-\j ^*°F(k)A° F (k + ») = -— (!>-- + -) - (20) 



3 Results for a massless theory with quartic interaction 

When J7(0) = A</> 4 /4!, the classical equation of motion is 

-d 2 T Mr) + ^l(r) = 0, 

C (O) = c (/3) = 0o ■ (21) 
The solution of this equation with the required boundary condition is given by 

Mr) = nc (ipt(r - /3/2), , (22) 

where nc is one of the twelve Jacobi Elliptic Functions [19], and ip t is defined implicitly by 
the following equation: 

^° = \[j nc(^/3/2, l/>/2) . (23) 
Substituting (22) in Eq. (15), we obtain 



di-k 1 



<p 2 t nc 2 (^(r-/V2),l/\/2)] ?7(r,fc 2 ) =0, (24a) 



r ? (0;fe 2 ) = l, g(/3;fc 2 ) = . (24b) 

We solve Eq. (24) numerically, and use t](t, k 2 ) in Eq. (17) in order to obtain the one-loop 
effective potential in the BET approach. 

The BET effective potential (in units of T 4 ) is plotted in Fig. 2 as a function of the 
dimensionless field P<f>o for a coupling constant A = 10 (the results for A = 1 are exhibited 
in Fig. 3). The result is displayed in the form of a band corresponding to a variation 
of the renormalization scale within the interval irT < [i < 4nT. One can see that the 
residual sensitivity to the renormalization scale is fairly small, suggesting that higher-order 
corrections are well under control. In Fig. 3, which shows the results for A = 1, this band is 
so narrow that it appears as a line. 

Our result is first compared to the classical potential itself, U(<fro)- One can see 
that it differs from the classical potential both at small field, due to the appearance of a 
quadratic mass term, and at large field due to large non-linear corrections. At small field, 
the standard one-loop result and BET are in good agreement, since both incorporate the 
effect of the thermal mass m 2 = AT 2 /24. However, for larger values of the field, differences 
appear. In fact, the one-loop result has an asymptotic behavior close to (A/24)(/3</> ) 4 (the 
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Figure 2: Plot of (3 4 V e g(fa) for a massless quartic interaction with A = 10. The shaded band 
represents the BET result, with a renormalization scale varying in the range [i € [7rT, 4ttT]. 
The thin solid line is the classical potential, i.e. U(fa) = X4>q/4\ in this case. The open 
circles represent the effective potential obtained from the classical action evaluated at the 
solution fa of the Euler-Lagrange equation with boundary value fa. The open triangles 
represent the standard 1-loop result, evaluated at fx = 2ttT. 



one-loop correction becomes very small at large field, and one simply recovers the classical 
potential), while the BET result approaches y^A/27 ((3 fa) 3 . One can see that this result is 
in fact dictated by the behavior of the classical solution fa[fa\: when the boundary field 
fa(x) — fa is large, the non-linear term in the classical Euler-Lagrange equation is very 
important, and the solution fa becomes a strongly non-linear function of the boundary 
field value fa. This non-linearity alters significantly the behavior of the BET effective 
potential at large fa. In fact, one can see that in this regime the BET effective potential 
is well approximated by the classical action evaluated at the solution fa. This suggests 
that for large fields, the main effect comes from the infinite sum of the non-linear tree level 
contributions. 



4 Conclusions 



The boundary effective theory (BET) framework provides a way to control the infrared 
divergences of thermal field theory in a well-defined and relatively simple way [7] . Previously, 
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Figure 3: Plot of V c g (<f>a) / j3 A for a massless quartic interaction with A = 1. Definitions are 
the same as in Fig. 2. 

we had computed the pressure of a massless hot scalar A<^ 4 theory [8], obtaining excellent 
agreement with up-to-date results from weak-coupling and screening perturbation theory 
[9, 10, 11]. In this paper, we have applied this method to the computation of the one- 
loop effective potential, following our previous work on the semiclassical thermodynamics of 
scalar fields [6]. 

The effective potential obtained within the BET formalism perfectly reproduces the 
standard one-loop result [2] for small fields, as expected, since the BET effective action 
contains very naturally the effect of the thermal mass. For large fields, BET goes beyond 
by incorporating the nonlinear corrections that become more and more important and are 
not captured by the standard one-loop calculation. We have also shown that our results are 
very stable with respect to variations of the renormalization scale, signaling a good behavior 
of the (semiclassical) series. 

A natural, and very useful, extension of this work would be treating the case of 
thermal symmetry restoration in the case of a double-well classical potential for the scalar 
field, with its consequences for spontaneous symmetry breaking and the description of phase 
transitions. However, this is not a straightforward extension, since the case of multiple wells 
presents nontrivial features related to the appearance of caustics and complex trajectories 
in the calculation of the semiclassical density matrix (see Refs. [20, 21] for a discussion). 
Nevertheless, we believe that the nonlinear corrections captured by the BET approach can 
be of great relevance in the description of phase transitions. Results in this direction will 
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be reported in the future [22] . 
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